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На примере РБМК-ЮОО произведено определение ресурса реакторного графита с учетом неравномерности энерговыделения по 
радиусу и высоте активной зоны, а также по радиусу графитового блока. Приведены: методика обработки массивов эксплуата- 
ционных параметров промышленного уран-графитового реактора; распределения объемной плотности тепловыделения с уче- 
том радиального и аксиального профилирования активной зоны ; распределения величины ресурса реакторного графита по ак- 
тивной зоне и графитовому блоку. 


Введение 

Анализ современного состояния исследований 
свойств реакторного графита и его работоспособ- 
ности в водографитовых реакторах позволяет 
сформулировать следующие выводы: 

1. Для корректного определения ресурса графита 
необходимо учитывать влияние сопутствующе- 
го у-излучения [1]. 

2. Плотность потока у-излучения зависит от места 
расположения графитового блока в активной 
зоне, от места в графитовом блоке и типа само- 
го блока. Поэтому необходимо проведение бо- 
лее точных расчетов полей у-излучения, тем бо- 
лее, что в каждой точке блока наблюдается раз- 
ная плотность потока быстрых нейтронов и эк- 
вивалентная температура облучения [2]. 


3. Для строгого определения ресурса графита 
необходимо учитывать, что нейтронно- и те- 
плофизические условия эксплуатации зависят 
не только от места расположения графита в ак- 
тивной зоне. Они являются немонотонными 
функциями времени и значительно изменяются 
в течение кампании. 

4. Использование средних интегральных значе- 
ний плотности потока быстрых нейтронов и 
температур при определении момента достиже- 
ния критического флюенса приводит к значи- 
тельной погрешности. 

В связи с этим в данной работе была поставлена 
задача: разработать метод определения ресурса ре- 
акторного графита в различных местах кладки, ис- 
ходя из параметров эксплуатации реального ядер- 
ного реактора. 


Технические науки 


Методика обработки массивов 

эксплуатационных параметров 

Одними из основных контролируемых параме- 
тров является мощность ядерного реактора в тече- 
ние определенного промежутка времени. Для 
упрощения работы с массивами эксплуатационных 
параметров необходимо провести их сверку. 

Исходными массивами являются значения те- 
пловой мощности реактора /V/ и длительности ми- 
крокампаний А//, где і - номер микрокампании, 
входящей в состав у'-ой кампании. Пример массивов 
эксплуатационных параметров приведен на рис. 1. 

По принципу прямого сравнения производятся 
выборка и затем свертка: 

Л?=Л?, где тФіг, 

Л^”=/Ѵ,\ где ЫІ, 
п,т= 1,2,...,/; к,1=\,2,...,І. 

Если где шп, то И' к п =^=..-И ш -Аі м = Д 

(к,т) 

Ѵ,п) 

если N"=N 1 , где Ы, то N”=7^"=...=^; Дг„, 2 = X Аі т> 

(к,т) 

(/.») 


и т.д. для \ѵ 3 ,ѵг 4 ,..., 

где /- количество микрокампаний; У - количество 
кампаний; - идентификатор микрокам- 

паний, в течение которых значения энергонапря- 
женности реактора совпадают. 


В результате свертки формируются два одно- 
мерных массива: И К и Аі„ где УѴ„, — тепловая мощ- 
ность реактора за промежуток времени Л? и „ 

Результаты показали, что погрешность опреде- 
ления флюенса по среднему значению тепловой 
мощности за год достигает 14%. 


Радиальное и аксиальное 
распределение температурных полей 


На основе анализа экспериментальных данных 
аппроксимированы радиальное и аксиальное ра- 
спределения объемной плотности тепловыделения 
с учетом радиального и аксиального профилирова- 
ния активной зоны (рис. 2). Коэффициент нерав- 
номерности энерговыделения по радиусу составля- 
ет кр 1,1, по колонне - к^ 1,3. Ширина плато со- 
ставила: по радиусу 3,7 м, по высоте 3,4 м. 

Максимальное и минимальное значения энерго- 
выделения в колонне определяются из соотношений: 


дтшах 

^ ' КОЛ 


Ы-к г 

9 

У1 

,1 ТК 


м тіп 

кол 


0,3 -К 


тах 
кол ? 


где N — средняя мощность ядерного реактора; 
к г - коэффициент неравномерности энерговыделе- 
ния по радиусу активной зоны с учетом радиального 
профилирования; п ТК - число топливных каналов. 

Расчетным путем установлено аксиальное ра- 
спределение температуры теплоносителя, при этом 
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максимальное значение объемной плотности тепло- 
выделения в топливе определялось из соотношения: 


% = ■ 


Коп ~ к - 


,-п ■ К тт ■ Н АЗ 


где к. - коэффициент неравномерности энерговы- 
деления по высоте активной зоны с учетом акси- 
ального профилирования; п шм - количество твэлов 
в топливном канале; К п - радиус топливной та- 
блетки; Н АЗ - высота активной зоны. 

Минимальное значение объемной плотности те- 
пловыделения на краях активной зоны составляет 
30 % от максимального значения и имеет место на рас- 
стоянии 0,5 м с каждой стороны по высоте колонны. 

Средняя температура теплоносителя на /-ом 
участке определяется соотношением: 

Т'тн = Т‘, и ' + <2 і /(С • с Р ), 


где / - номер участка по высоте колонны, /= 0 соот- 
ветствует входу теплоносителя в колонну; 
0=д‘- л- Яр Т -Аіп тюі , - количество тепла, выделяемое 
на /-ом участке, где - среднее объемное тепловы- 
деление в топливе на /-ом участке; А г - длина 
участка; С - расход теплоносителя; с~ р - удельная 
теплоемкость теплоносителя, п шм - количество 
твэлов в технологическом канале. 


Определено, что объемная плотность тепловы- 
деления в графитовых блоках прямо пропорцио- 
нальна объемной плотности тепловыделения в то- 
пливе и определяется соотношением: 


Ч, 


_ ОР'зХр 8 Т0П д і 
Е / ^ ' Е ЛЗ ' 


где Е„ - среднее значение энергии нейтронов деле- 
ния; Е / — энергия деления II 235 ; (^Е 5 ) гр , замед- 
ляющие способности графита и активной зоны, со- 
ответственно; Е топ - площадь поперечного сечения 
(вид сверху), занимаемая топливом в активной зо- 
не; - площадь поперечного сечения активной 
зоны; Д=4...6 - характеризует отношение г/, (за счет 
у- квантов и нейтронов)/#,, (за счет нейтронов). 


Температура на внутренней поверхности графи- 
тового блока на /-ом участке определялась из соот- 
ношения: 


Т^ш)~ Т-іи +ЧІ 


с 

°тр р 
+ -у + К за 

'"тр 


Л 

2 


где Ф 


2 -л ■ К ш 


- поверхностная плотность те- 


плового потока на /-ом участке; здесь 5 гр - площадь 
поперечного сечения ячейки, занятая графитом; 
К тк - радиус технологического канала; а - коэф- 
фициент теплоотдачи; 5 тр - толщина трубы техно- 
логического канала; Х щ - коэффициент теплопро- 
водности для материала трубы технологического 
канала; К ш - термическое сопротивление зазора, 
заполненного азотно-гелиевой смесью. 


Графит, подвергаясь радиационному поврежде- 
нию, значительно деформируется, в результате че- 
го в одних местах величина газового зазора умень- 
шается, в других - увеличивается. В силу того, что 
коэффициент теплопроводности газа очень мал, 
небольшие изменения зазора приводят к значи- 
тельным изменениям термического сопротивле- 
ния. Анализ экспериментальных данных показал, 
что среднее значение термического сопротивления 
газового зазора составляет 200/^ шах , где #, тах - 
максимальное значение поверхностной плотности 
теплового потока. 


Средняя эквивалентная температура по сече- 
нию графитового блока на /-ом участке определя- 
ется из соотношения: 


Ті = Т,(К тк ) + 0,Ш5 



(А -АХ 


где Я - коэффициент теплопроводности графита; 

А =1і/ Кэ 1-й = 0,25 ^ э -Ак где к _ эквива- 
1 т {2К тк )’ Л і 2 К 2 Э ’ 

лентный_ геоме- 
трический радиус блока графита; К э =а/^к, где а - 
сторона графитового блока. 



Рис. 2. Радиальное и аксиальное распределения объемной плотности тепловыделения (темная линия - результат линейной 
аппроксимации) 


Технические науки 


Определение ресурса реакторного графита 

Значение критического флюенса является силь- 
ной функцией эквивалентной температуры облу- 
чения. Так, в условиях РБМК-1000, функциональ- 
ная зависимость критического флюенса удовлетво- 
рительно описывается полиномом 2-ой степени: 

Р 0 (Т) = а 0 + а х -Т + а 2 - Т 2 , 

где я 0 = 5,92; а,=-9,0610 _3 ; д 2 =3,ЗН0“ 6 - при темпе- 
ратуре более 300 °С; а 0 =О,322; а { =— 1,9 МО -3 ; 
й 2 =1,0210“ 5 - при температуре до 280 °С; Т измеря- 
ется в градусах Цельсия; Т 0 - в ІО 22 см 2 . Значения 
критического флюенса при температурах от 280 до 
300 °С получены путем линейной аппроксимации. 
При определении коэффициентов полинома были 
использованы экспериментальные пары значений 
Т 0 и Г для графита ГР-280 при плотности потока со- 
путствующего у- излучения Ф =5Т0 14 см 2 с~‘ [2]. 

Согласно расчетам критический флюенс на са- 
мом энергонапряженном участке для внутренней 
поверхности графитового блока составил 
2,3' ІО 22 см -2 , его среднее значение по графитовому 
блоку составило 1,65- ІО 22 см -2 . 

Значения плотности потока быстрых нейтронов 
на каждом участке получены из соотношения: 

со 

ч:=М яд -А- | а ы [Е)- 8ы (Е)-Е-Ф 1 [Е)0Е = 

ШкэВ 

= Ы яд -А -аеі-ё еІ -Ё-Фі, 

где Ы яд - концентрация ядер графита; о;, - микро- 
скопическое сечение упругого рассеяния; § е1 - доля 
энергии, теряемая нейтроном при упругом рассея- 
нии; Е — энергия нейтронов; Ф, - плотность потока 
быстрых нейтронов на /-ом участке; а л Ф,— значе- 
ния. полученные при средней энергии нейтронов Е . 

Результат расчета для РБМК-1000 при работе на 
номинальной мощности 272 эффективных суток в 


Таблица 1. Распределение ресурса реакторного графита (в 
годаах) по активной зоне; г - аксиальная; г ~ ра- 
диальная координаты 


Г 

і X. 
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140 

140 

140 

140 

140 
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233 
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114 

114 

114 

114 
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191 

233 

410 

3 

/1 

71 

/1 

/1 

71 

8/ 
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26/ 

4 

44 

44 

44 

44 

44 

55 

78 

97 
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5 

29 

29 

29 

29 

29 

37 

53 

67 

123 

6 

23 

23 

23 

23 

23 

29 

43 

55 

102 

7 

22 

22 

22 

22 

22 

28 

42 

53 

100 

8 

22 

22 

22 

22 

22 

28 

42 

53 

100 

9 

22 

22 

22 

22 

22 

28 

42 

53 

99 

10 

22 

22 

22 

22 

22 

28 

42 

53 

99 

11 

21 

21 

21 

21 

21 

28 

41 

52 

99 

12 

21 

21 

21 

21 

21 

28 

41 

52 

99 

13 

21 

21 

21 

21 

21 

2/ 

41 

52 

99 

14 

21 

21 

21 

21 

21 

27 

41 

52 

98 

15 

21 

21 

21 

21 

21 

28 

42 

53 

100 

16 

2/ 

2/ 

2/ 

2/ 

2/ 

34 

51 

64 

120 

17 

41 

41 

41 

41 

41 

51 

74 

93 

171 

18 

66 

66 

66 

66 

66 

82 

116 

144 

261 

19 

105 

105 

105 

105 

105 

129 

182 

224 

401 

20 

134 

134 

134 

134 

134 

159 

217 

275 

467 


течение года дан в табл. 1. В таблице приведен ре- 
сурс реакторного графита (по радиусу активная зо- 
на разбита на 9 участков, нулевой участок соответ- 
ствует центру активной зоны; по высоте — на 
20 участков, первый - входу теплоносителя в ак- 
тивную зону). При использовании среднего инте- 
грального значения температуры ресурс для всего 
реакторного графита составил около 50 лет. 

Изменение потока повреждающих нейтронов 

по радиусу ячейки 

Эквивалентная температура облучения графита 
(выше 300 °С) возрастает по радиусу ячейки, следо- 
вательно, значение критического флюенса снижа- 
ется. При этом также снижается значение плотно- 
сти потока повреждающих нейтронов. Имеет 
смысл рассмотреть влияние неравномерности 
плотности потока быстрых нейтронов по ячейке на 
ресурс. 

Эквивалентная температура в графитовом бло- 
ке определяется соотношением: 

Ті(г) = Т і (Я тк ) + 

I _ <- К тк У г г 2 -Вф к Л 
Я-(1 -К 2 тк /К 2 э ){ к тк 2- К 2 У 

где г - радиальная координата, К тк <г<К э . 



Рис. 3. Изменение плотности потока быстрых нейтронов по 
радиусу ячейки 


Величина плотности потока быстрых нейтро- 
нов по сечению технологического канала изменя- 
ется слабо, поэтому ее можно считать постоянной. 
Выбранное распределение представлено на рис. 3 и 
описывается соотношением: 

Ф й =Фо-ехр[-Е,-(г-і^)], 

где Ф 0 - значение плотности потока быстрых ней- 
тронов в сечении технологического канала, опре- 
деляемое из известного среднего по ячейке значе- 
ния; X, - полное макроскопическое сечение для 
энергии нейтронов 2,5 МэВ. 
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Результаты расчета радиуса для центральной 
ячейки реактора РБМК-1000 приведены в табл. 2. 


Таблица 2. Распределение ресурса реакторного графита по 
радиусу эквивалентной ячейки 


Радиус, см 

Ресурс, лет 

4,40 

20,8 

5,36 

21,6 

6,32 

22,9 

7,28 

24,5 

8,24 

26,5 

9,20 

28,9 

10,16 

31,9 

11,12 

35,4 

12,08 

39,6 

13,04 

44,6 

14,00 

50,5 


Заключение 

В экспериментальных оценках ресурса реактор- 
ного графита распределение энерговыделения по 
радиусу и высоте активной зоны является средней 
интегральной характеристикой. Значение ресурса 
при таком подходе для уран-графитовых реакторов 
определено в 35 лет. Эти оценки носят противоре- 
чивый характер. Например, графитовая кладка про- 
мышленных уран-графитовых реакторов, имеющих 
больший «возраст», оказывается менее поврежден- 
ной, чем в реакторах, введенных в эксплуатацию 
позже. При ресурсе 35 лет реакторы продолжают ус- 
пешно функционировать в течение 50 лет, а при 
тщательном исследовании графитовой кладки ока- 
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зывается, что она может эксплуатироваться еще 
около 10 лет. Таким образом, расчет ресурса необхо- 
димо проводить для каждой отдельной колонны с 
учетом её местоположения по отношению к орга- 
нам управления и защиты. Отличие расчетных зна- 
чений эквивалентной температуры облучения и ре- 
сурса реакторного графита в различных точках ак- 
тивной зоны от усредненных может достигать 23 %. 

Оценки ресурса, как правило, проводятся с ис- 
пользованием экспериментальных зависимостей 
критического флюенса от температуры облучения, 
полученных в исследовательских реакторах. При 
этом под температурой облучения не всегда подра- 
зумевается эквивалентная температура, обусло- 
вленная взаимодействием нейтронов и у-квантов с 
графитом. Результаты расчетов показали, что по- 
грешность определения критического флюенса, 
вызванная данным обстоятельством и использова- 
нием среднего значения энергонапряженности 
кладки, за год может достигать 14 %. 

С целью увеличения степени адекватности рас- 
четных данных ожидаемым значениям следует учи- 
тывать влияние на величину критического флюен- 
са действующего значения плотности потока пов- 
реждающих нейтронов и его спектра. 

Использование усредненных во времени и по 
активной зоне микро- и макроскопических пара- 
метров эксплуатации графита во многих случаях 
приводит к занижению реального ресурса. 

2. Карпухин В.И., Николаенко В.А., Кузнецов В.Н. Критический 
флюенс нейтронов как фактор, определяющий ресурс графитовой 
кладки РБМК //Атомная энергия. - 1997. - № 5. - С. 325-330. 
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